Fluid adsorption near an apex: Covariance between complete and critical wetting 
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Critical wetting is an elusive phenomenon for solid-fluid interfaces. Using interfacial models we 
show that the diverging length scales, which characterize complete wetting at an apex, precisely 
mimic critical wetting with the apex angle behaving as the contact angle. Transfer matrix, renor- 
malization group (RG) and mean field analysis (MF) shows this covariance is obeyed in 2D, 3D and 
for long and short ranged forces. This connection should be experimentally accesible and provides 
a means of checking theoretical predictions for critical wetting. 

PACS numbers: 68.08.Bc, 05.70.Np, 47.20.-k, 68.35.Md 
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Advances in the controlled fabrication of micropat- 
terned substrates have stimulated the experimental and 
theoretical study of fluid adsorption at tailored surfaces 
[0, H ||, ^, ||. For example, Mistura and co-workers ||] 
have recently investigated complete wetting of Ar on sev- 
eral parallel arrays of wedges and apexes. They show con- 
vincingly that the adsorption within the (independent) 
wedge regions is geometry dominated and distinct from 
the planar complete wetting case. As well as having im- 
plications for microfluidics such studies have also revealed 
a number of unexpected results relating interfacial fluc- 
tuation effects and substrate geometry which have wide 
application to other phase transitions (see later). Here 
we use effective Hamiltonian theory to show that com- 
plete wetting on apex shaped substrates reveals a hidden 
connection (covariance) with critical (continuous) wet- 
ting transition occurring on planar surfaces ^, 0|. The 
covariance emerges when one considers how, at bulk co- 
existence, the mean height ^^(ck) of the unbinding inter- 
face above the apex tip depends on the apex angle a. 
For shallow apexes we show that Ia{o) is identical to the 
mean interfacial thickness occurring at a particular class 
of critical wetting transition with the apex angle playing 
the role of an effective contact angle. The covariance is 
valid for 2D and 3D apexes, for arbitrary intermolecular 
forces and is, we believe, a general feature which should 
also be present in more microscopic models. 

The central result of our article is the following covari- 
ance relation for the interfacial height probability dis- 
tribution function (PDF) which precisely quantifies the 
influence of the apex geometry on the complete wetting 
film. We emphasize that the PDF contains a great deal 
of information determining the (local) interfacial height, 
roughness and the scaling properties of the density pro- 
file. Let P^{1, a) denote the PDF for the interface height 
above the apex where the superscript refers to the repul- 
sive binding potential U{1) ~ 31"^ for a complete wet- 
ting film (see below). Let [1; 9) denote the PDF for a 
planar critical wetting transition written in terms of the 
contact angle 9. Here V{1) is the binding potential for 
a critical wetting transition which, of necessity, contains 
attractive and repulsive terms. The covariance relation. 
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valid for small a and at bulk coexistence, reads 

P^{l,a)^P^{l-a) 

where the covariant binding potential 

V{1) ^ -Al-^'"^ + U(l) 

with X = niin(p, r) and A determining the effective con- 
tact angle (related to a). Here r = 2(1 — C)/C and C, 
are the entropic repulsion and interfacial wandering ex- 
ponents respectively . Thus the apex locally binds the 
complete wetting film and induces an effective attractive 
term in the binding potential, twice the range of the dom- 
inant intermolecular or entropic repulsive term. A similar 
rule applies in 3D for exponentially decaying potentials. 
Recall that in contrast to abundant experimental studies 
of complete wetting (including recent work on systems 
with short ranged forces [^) critical wetting transitions 
are rather rare |^ for which no examples are known 
currently for solid- ffuid interfaces. The covariance dis- 
cussed here provides a means of effectively inducing crit- 
ical wetting behaviour using complete wetting films. 

Consider the interface between an infinite apex and 
a bulk vapor at temperature T and chemical potential 
A* < fJ-satiT) (see Fig. |l|). We suppose that the flat 
wall (a = 0) is completely wet by the liquid phase at 
coexistence J/i = fisat ~ l^{T) = 0"*" corresponding to zero 
contact angle. The wall shape is described by a height 
function za = — tana|a;| in the (x, z) plane although we 
shall only be interested in the case of shallow apexes for 
which we may approximate tana ~ a. Macroscopically 
far from the apex tip the height of the interface above the 
wall is the same as that occurring for a flat wall. Since 
the liquid-vapor interface is required to round the apex, 
surface tension restrictions imply that the local height I a 
above the apex tip is smaller and remains flnite even in 
the limit of bulk coexistence. We wish to evaluate the 
mean interfacial height Za(ck) and interfacial roughness 
(r.m.s. interfacial width) S,a{cx) at bulk coexistence and 
the critical exponents 

-Pa 
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For the 3D apex we also wish to determine the transverse 
correlation length ^y(a) ^ a~'^^ pertinent to correlations 
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FIG. 1: Schematic illustration of 3D apex complete wetting, 
showing a section of a typical interfacial configuration above 
the tip. Diverging length scales are highlighted. 



along the apex. Correlations in the x direction are not 
described by a finite correlation length and fluctuations 
are not localised to a region near the apex. 

We begin with the 2D apex. The starting point for our 
calculations is the interfacial Hamiltonian model 



U{l + a\x\) 



(4) 



where l{x) is the local height of the interface above the 
z — Q reference line, S is the (reduced) stiffness coeffi- 
cient (surface tension) of the liquid-vapor interface and 
U{1) denotes the binding potential modelling the com- 
plete wetting behavior pertinent to the planar system 
a = 0. For shallow apexes, it is permissable to assume 
the interface interaction with the wall occurs via the rel- 
ative vertical height I = / + a|a;|. The binding potential 
U{1) has an infinite hard wall repulsion and decays as 



[/(/) = [pi - Pg)l3S^il + Bl-P 



(5) 



where, for the moment we have allowed for a finite bulk- 
order field > 0. Here i3 is a positive Hamaker con- 
stant whilst p accounts for the range of the intermolec- 
ular forces. Exponentially decaying binding potentials 
will also be considered in our discussion of the 3D apex. 
For dimensions in which the free liquid- vapor interface is 
rough, critical effects at planar complete wetting transi- 
tions fall into two classes |l^ : a MF regime for p < r 
and a fluctuation-dominated regime for p > t. Heuristi- 
cally this arises from the interplay between the direct in- 
temolecular repulsion ^ 1~p and the effective entropic re- 
plusion U fi{l) ^ In the present paper we restrict our 
attention to pure systems for which the interface is rough 
for d < 3 and C, = (3 — (i)/2. For fixed p the upper critical 
dimension for complete wetting is dco = 3 — 4/ (p+ 2) [ p^ . 

The model can be studied using transfer matrix meth- 
ods previously developed for the wedge geometry [pT|. 



Care must be taken in defining an infinite apex geometry 
and it is convenient to first consider a finite apex ex- 
tending over the range [— L/2,L/2] and impose periodic 
boundary conditions at the edges. After taking the ther- 
modynamic limit L — > c» at finite 5p > Q it is straight- 
forward to derive an expression for the interfacial height 
probability distribution function at arbitrary position x 
along the wall. At the apex mid-point symmetry consid- 
erations simplify this expression considerably 



P^(/;a)(x|^o(OI' 



(6) 



where denotes the ground-state wave function solv- 
ing the Schrodinger equation 
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with boundary conditions ip{Q) = ipi^o) — 0. We now 
focus on the complete wetting limit 6p, 0. Macro- 
scopically far from the apex the interface unbinds from 
the wall. Close to the apex however the interface re- 
mains bound due to the pinning exponential term in (|^). 
As a — + 0, three different critical behaviors are found: 
(I) A MF regime for p < 2, characterized by Gaussian 
fluctuations with Ia{(^) ^ ^a(ck); (II) a marginal case 
for p — 2; and (III) a fluctuation dominated regime for 
p > 2 with universal critical behavior and large scale 
fluctuations /a (a) ~ ^A{ct) ~ a~^. The explicit expres- 
sions for the large distance/scaling behavior of P^{1; a), 
determining the critical singularities, are given by: 



exp (-2M 



4\/2SB ;l-p/2\ 
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PX{l;a)^ { /i+v^+8^exp(-2SaO 



Pexp(-2Ea/) 



p<2 



p>2 



(8) 

In the MF regime, a saddle point evaluation reveals that: 



lA{a) 



f 2B V'" 



(9) 



and va — 1/2. The critical exponents are con- 

tinuous at p = 2. These results completely classify the 
asymptotic critical behavior for complete wetting at a 
2D apex in pure systems. At this point we make two 
remarks: 

A: The values of the critical exponents follow from a 
simple mean-fleld/entropic repulsion argument. Ignoring 
fluctuations the equilibrium interfacial proflle is obtained 
from minimization of the effective Hamiltonian. A first 
integral of the Euler-Lagrange equation determines the 
mid-point height at bulk coexistence according to 



U{Ia) 



(10) 
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and leads directly to the result (||) valid in the MF regime 
(p < t). For p > T interfacial wandering leads to an 
entropic repulsion Ufi ^ l^^ . Thus we should expect two 
regimes with (3a = max(2/p, C/(1~C)) in agreement with 
the exphcit calculation for C = 1/2. For later purposes 
observe that the MF equation is also appropriate for 
higher dimensional apexes. 

B: The PDF's and associated critical exponents are 
identical to those occurring at a certain class of 2D criti- 
cal wetting transition. At a critical wetting transition the 
mean height of the interface Z^, roughness and par- 
allel correlation length for a planar substrate diverge 
as the temperature (say) is increased towards a wetting 
temperature (/x = fJ-sat)- This is equivalent to the 
contact angle of a sessile drop vanishing as T ^ . 
The standard interfacial model for this is 

WZ]^/d.||(|)%no} (11) 

where V{1) denotes an appropriate binding potential. 
The associated PDF P^{l;e) can be calculated using 
standard methods which map the problem onto one di- 
mensional quantum mechanics [lO| . In particular con- 
sider 2D critical wetting transitions occurring for the 
class of potentials (||) with x = min(p, 2). For such po- 
tentials calculations show that A a 9. A straightforward 
calculation of the PDF's P^{1\9) for p < 2, p = 2 and 
p > 2 yields results that, in the critical limit (small 9) arc 
identical to (^ provided we set 9 — a |jl^. Apex com- 
plete wetting precisely mimics the properties of a critical 
wetting transition. From the covariance of the PDF's it 
follows that the mean-interfacial heights satisfy 

lA{a) = lZ{a) (12) 

where the RHS is understood to represent the mean 
height at a planar critical wetting transition with the 
covariant potential (H). It is notable that apex covari- 
ance is obeyed at MF level and beyond, and therefore not 
necessarily related to fluctuation induced effects such as 
hyperscaling. In particular the relation ( [l^ ) follows di- 
rectly from comparing the solution of the MF equation 
( p^ ) with the position of the minimum of the critical wet- 
ting potential (||) (with p). The covariance relations 
for apex complete wetting are similar, but not identi- 
cal, to those which exists for 2D wedge filling transitions. 
For both pure and impure systems 2D filling transitions 
mimic the properties of planar critical wetting transitions 
with short-ranged forces in contrast to the present apex 
problem where the equivalent critical wetting transition 
has long(er)-ranged forces. Wedge covariance and fill- 
ing are closely related to the Indekeu-Robledo conjecture 
for the line tension |l3[ , and the unzipping transition for 
stranded polymer chains p^ . It is likely that similar 
connections may also apply for the apex geometry. 



Having discussed the 2D apex in detail it is straight- 
forward to generalise the results to 3D systems based on 
the interfacial model 

I3Ha[1]^ y"dx||(VO' + [/(? + a|x|)| (13) 

where for purposes of generality we have written x ~ 
(x, X|| ) where xy denotes the d — 2 dimensional vector 
along the apex axis. In 3D x = {x, y) but it is instructive 
to also consider the generalized apex for 2 < d < 3 since 
this gives clear indication that the covariance relations 
extend to higher dimensions. First rewrite the Hamilto- 
nian in terms of the relative height I. The critical behav- 
ior follows from elementary RG considerations. Under 
rescaling x ^ x' = x/6, = Ih^'^ the the renor- 

malised tilt angle and Hamaker constant are a' = ah^^'^ 
and B' = B\P'~^p~'^'^ respectively. Thus a is always a 
relevant scaling field whilst the intermolecular forces are 
only relevant for p < t. The criticality falls into two 
scaling regimes consistent with the explicit 2D results: 

(I) A MF regime for p < r with [3 a = 2/p, Vy = 2/p-|- 1 
and VA = Qvy for which ( p^ is valid. 

(II) A fiuctuation-regime for p > t describing the uni- 
versality class of systems with short-ranged forces with 
Pa = C/(l - C), - 1/(1 - C) and va = C^y. 

Note that for fixed p the upper critical dimension dA — 
3 — 4/(p-f 2) and is unchanged from the planar complete 
wetting result dco- Remarks A and B made earlier about 
the 2D results also apply in higher dimensions, where the 
planar covariant effective Hamiltonian is: 

(3H^[l]^ J rfx||(V0' + V(0} (14) 

For the physically relevant case d — 3, MF theory 
is valid for all long-ranged intermolecular forces (finite 
p). Thus for non- retarded van der Waals forces we pre- 
dict ^(a) ~ y^2B/J^a^ and observe that this is identi- 
cal (covariant) with the growth of the interfacial thick- 
ness at a critical wetting transition with binding po- 
tential V — ~A/l -\- B/P. Similar remarks apply for 
the interfacial roughnesses at the respective transitions 
(^^(a) ^ ^± ^ %/— In a). Our final task is to address the 
issue of covariance for the marginal case of 3D systems 
with short-ranged forces. For this we use the interfacial 
model ( p^ with binding potential U — B exp{— nl) and k 
the inverse bulk correlation length. At MF level we find 
KlA{a) ^ — 21na whilst solution of the Ornstein-Zernike 
equation for the height-height correlation function along 
the apex tip yields ~ ||l^. Thus the MF expo- 
nents for short-range forces are f3A = 0(ln), Vy — 1 and 
are consistent with the ( ^ limit of the short-ranged 
exponents detailed in (II) . Beyond MF we anticipate that 
Vy is unchanged but that the logarithmic divergence of 
Ia is altered. To study this we employ the same ap- 
proximate linear functional RG approach used for the 3D 
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planar wetting transition but with an appropriately mod- 
ified matching condition for the new geometry [|6|. Un- 
der the rescaling of / and x described earlier (with C = 0) 
the binding potential maps as U{1) — *■ i?('''[J7] which ef- 
fectively coarse-grains the potential over the interfacial 
roughness. Since the angle and correlation length rescale 
as a' — ba and Cy — ^y/b one can curtail the renormal- 
ization at 6 ~ at which scale both the renormalized 
angle and transverse correlation length are of order unity 
and fluctuation effects are negligible. Matching with MF 
theory implies S/2 = R'--^^"'>[U{Ia)] and a simple calcu- 
lation yields 



KlA{a) ^ —(2 + Lu)liia 



(15) 



where uj = K^/(47rI]) is the usual wetting parameter and 
we have assumed w < 2 as pertinent to the bulk Ising 
universality class. This critical behavior is again consis- 
tent with covariance as can be seen by comparison with 
the 3D short-ranged critical wetting transition described 
by the model (14) with potential 



V{1) = -Ae 



Be 



■J>0 



(16) 



where, guided by our earlier findings, we have included 
an attractive term which is twice the range of the direct 
repulsion. This is the same binding potential appearing 
in the standard theory of short-ranged critical wetting 
except for a trivial factor of two in the definition of the 
inverse bulk correlation length. Accordingly, rescaling k 
and u! in the known RG results for critical wetting |16| 



(17) 



behaviour of the critical wetting transverse correlation 
length written in terms of the contact angle S.\\{0) ''^ 9"^ . 
Interestingly one still finds the critical behavior ( |l5| ) 
if one improves the apex calculation to account for a 
position-dependent stiffness coefficient [|l^. That is, even 
if planar short-ranged wetting transitions are driven first- 
order by a stiffness instability mechanism, the apex still 
mimics the properties of critical wetting. 

We finish with comments relevant to experimental 
studies on periodic systems. Close to bulk coexistence 
the interfacial height above a single, infinite apex shows 
scaling behavior Ia ^ a 
ponent ~ 2 + (Sa 

wavelength L) finite size effects modifies this to Ia 



^'*W^(i)((5^a~^'*) with gap ex- 
On a periodic array (with 



W^(2)(<5Aia-^^,L/e^°), where ^ff 



is the 



Thus the adsorption above the apex tip only behaves 
like a single apex for sufficiently large L ^ ^y°. This is 
equivalent to the requirement that the vertical distance 
between the apex tip and the wedge trough, aL/2, is 
much larger that the local height of the interface above 
the wedge bottom 1^ ~ Ea"^ /2f}5fi{pi — Pg). 

In this paper we have shown that complete wetting at 
an apex mimics precisely planar critical wetting. Taken 
together with similar covariance relations for wedge fill- 
ing, there is clear evidence of a fundamental connection 
between contact and geometric angles. Further work is 
required to understand such covariances at a deeper level. 
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where the (rescaled) wetting parameter w < 2. Note 

also that the divergence of ~ is similar to the [7] 



[8] 

[9] 
[10] 

[11] 

[12] 

[13] 

[14] 

[15] 
[16] 
[17] 



transverse correlation length for planar complete wetting. 
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